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AGING STEM CELLS 


A Werner syndrome stem cell model 
unveils heterochromatin alterations 
as a driver of human aging 


Weiqi Zhang,* Jingyi Li,?* Keiichiro Suzuki,?* Jing Qu,** Ping Wang," Junzhi Zhou," 
Xiaomeng Liu,” Ruotong Ren,’ Xiuling Xu,’ Alejandro Ocampo,” Tingting Yuan,* 

Jiping Yang,’ Ying Li,’ Liang Shi,” Dee Guan,’ Huize Pan,’ Shunlei Duan,’ Zhichao Ding,’ 
Mo Li,” Fei Yi, Ruijun Bai,* Yayu Wang,” Chang Chen,’ Fuquan Yang,’ Xiaoyu Li,” 
Zimei Wang,® Emi Aizawa,” April Goebl,”? Rupa Devi Soligalla,® Pradeep Reddy,” 
Concepcion Rodriguez Esteban,” Fuchou Tang,” ®™"?+ 

Guang-Hui Liu,?*™+ Juan Carlos Izpisua Belmonte*+{ 


Werner syndrome (WS) is a premature aging disorder caused by WRN protein deficiency. 
Here, we report on the generation of a human WS model in human embryonic stem 
cells (ESCs). Differentiation of WRN-null ESCs to mesenchymal stem cells (MSCs) 
recapitulates features of premature cellular aging, a global loss of H3K9me3, and 
changes in heterochromatin architecture. We show that WRN associates with 
heterochromatin proteins SUV39H1 and HPle and nuclear lamina—heterochromatin 
anchoring protein LAP2f. Targeted knock-in of catalytically inactive SUV39H1 in wild-type 
MSCs recapitulates accelerated cellular senescence, resembling WRN-deficient MSCs. 
Moreover, decrease in WRN and heterochromatin marks are detected in MSCs from 
older individuals. Our observations uncover a role for WRN in maintaining heterochromatin 
stability and highlight heterochromatin disorganization as a potential determinant of 


human aging. 


erner syndrome (WS), also known as 
adult progeria, recapitulates certain as- 
pects of human physiological aging 
(1). WS is caused by mutations in the 
WRN gene, resulting in loss of WRN 
expression or function (J). WRN protein plays 
roles in DNA replication, transcription, repair, 
and recombination, as well as telomere main- 
tenance, indicating that one of the major causes 
for WS pathogenesis relates to genomic insta- 
bility (7, 2). Epigenetic alterations have been 
associated with cellular aging in diverse model 
organisms (2-4). In humans, somatic cells derived 
from patients with premature aging syndromes 
are characterized by loss of heterochromatin 
marks (5-7). However, it is unclear whether 
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epigenetic dysregulation is involved in WS 
pathogenesis. 

Generation of patient-specific induced pluri- 
potent stem cells (iPSCs) represents a promising 
avenue to model and study human aging and 
aging-associated disorders (8). WS-specific iPSC 
lines may constitute an ideal source for in vitro 
modeling of WS. However, we found that WS 
patient fibroblast lines deposited in different cell 
banks presented severe karyotypic abnormalities 
and secondary DNA mutations associated with 
advanced stages of WS pathology. To create an 
unbiased human WS cellular model, we sought 
to generate an isogenic WS embryonic stem cell 
(ESC) line by knocking out exons 15 and 16 of the 
WRN gene encoding the conserved DNA helicase 


domain (9). After two rounds of homologous re- 
combination using helper-dependent adenoviral 
vector (HDAdV) (10, 11), we successfully generated 
homozygous WRN-null ESC lines (ESCs-WRN “) 
(Fig. 1, A and B, and fig. S1, A to D). ESCs-WRN~/- 
expressed pluripotency markers, maintained nor- 
mal karyotype, and could differentiate into all 
three germ layers (Fig. 1A and fig. S2, A to E). 
ESCs-WRN~’~ lacked detectable WRN protein, as 
determined by Western blot using antibodies spe- 
cific to the N or C terminus of WRN (Fig. 1B). No 
difference in cell cycle kinetics and cell growth 
rate between wild-type and WRN-null ESCs was 
observed (fig. S2, F to H). 

WS patients are mainly characterized by pre- 
mature aging pathologies associated with degen- 
eration of mesodermal tissues, i.e., osteoporosis, 
atherosclerosis, and gray hair (7). We hypothesized 
that WS patients may suffer from an accelerated 
exhaustion of the mesenchymal stem cell (MSC) 
pool. This was tested by differentiating ESCs- 
WRN ^ into MSCs. MSCs-WRN ~~ expressed MSC- 
specific cell surface markers CD73, CD90, CD105; 
lacked expression of MSC-irrelevant antigens, in- 
cluding CD45, CD34, and CD43 (fig. S3A); and 
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could differentiate into osteoblasts, chondrocytes, 
and adipocytes (fig. S3, B and C) (12). 

Upon serial passaging, WRN-deficient MSCs re- 
capitulated major phenotypes of premature aging, 
including premature loss of proliferative poten- 
tial, increased number of senescence-associated-B- 
galactosidase (SA-B-gal)-positive cells, up-regulated 
expression of aging-associated genes p16™*** and 
p21" and activation of senescence-associated 
secretory phenotype (SASP) (Fig. 1, C to E, and fig. 
S3, D to G) (13). Moreover, when WRN-deficient 
MSCs expressing luciferase were transplanted 
into the muscle of nonobese diabetic/severe com- 
bined immunodeficiency (NOD/SCID) mice, they 
underwent an accelerated attrition compared 
to wild-type MSCs (Fig. 1F and fig. S3H). These 
results demonstrated that the loss of WRN pro- 
motes premature senescence in MSCs. 

WRN deficiency in MSCs resulted in elevated 
DNA damage response (DDR), indicated by in- 
creased nuclear foci for 53BP1, y-H2AX, and phos- 
phorylated ATM/ATR substrates (fig. S4, A to C). 
Restoration of WRN activity by lentivirus-mediated 
expression in MSCs-WRN~"~ resulted in partial al- 
leviation of DDR and cellular senescence (fig. S4, 
D and E). To investigate potential chromosomal 
abnormalities resulting from the loss of WRN pro- 
tein, we performed genome-wide copy number 
variation (CNV) analysis by deep sequencing. In 
the time frame examined, genomic integrity was 
minimally affected in MSCs-WRN~’~ (fig. S4F). 

Epigenetic alteration has been considered as a 
hallmark of aging (2). MSCs-WRN~/~ showed a 
distinct nuclear Hoechst 33342 staining pattern 
with markedly enlarged nuclei and a high pixel-to- 
pixel coefficient of variation (CV) value, indicating 
possible changes in chromatin structure (Fig. 2A 
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and fig. S5A). Moreover, WRN-deficient MSCs exhib- 
ited accelerated diminishment of heterochromatin- 
associated inner nuclear membrane (INM) proteins 
LAP2ß and LBR and reduced heterochromatin 
structure underneath the nuclear envelope, as in- 
dicated by immunostaining and electron micros- 
copy (Fig. 2B and fig. S5, B and C) (14). These 
results suggest a progressive disorganization of 
heterochromatin in WRN-deficient MSCs. 

Further investigation of heterochromatin re- 
organization at histone and DNA levels revealed 
marked down-regulation of the constitutive het- 
erochromatin mark H3K9me3 (trimethylated 
histone H3 at lysine-9) in MSC-WRN~‘~ (Fig. 2C 
and fig. S5, D and E). In contrast, H3K27me3 showed 
slight down-regulation, whereas H3K4me3, a mark 
for euchromatin fiber, exhibited comparable levels 
between WRN-deficient and wild-type MSCs (Fig. 
2C and fig. S5, D and E). We did not observe ob- 
vious genome-wide alteration of 5-methylcytosine 
(5mC) in WRN-deficient MSCs (Fig. 2C). Bio- 
informatic analysis identified 73 H3K9me3- 
enriched “mountains” throughout the genome 
in MSCs-WRN*’*, which are characterized by 
>20 kb of consecutive peaks of H3K9me3 (Fig. 
2D). Of these H3K9me3 mountains, 28 (38%) 
were lost in MSCs-WRN ~~ (Fig. 2D). Interestingly, 
24 (86%) of these impaired H3K9me3 mountains 
resided in subtelomeric or subcentromeric re- 
gions (Fig. 2, D and E, and table S1). 

RNA sequencing (RNA-seq) identified 1047 
RefSeq genes that showed differential expression 
in MSCs-WRN~^ (fig. S6, A and B, and table S2). 
The most obviously down-regulated genes were 
centromere-packaging proteins and components 
of the nuclear membrane (fig. S6, A to E, and 
table S3). These results indicate alterations in 
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nuclear structure and epigenomic organization, 
potentially leading to a progressive loss of het- 
erochromatin structure in MSCs as a consequence 
of WRN depletion. 

In agreement with previous reports describing 
WRN as a telomere-associated protein required 
for telomere maintenance (15), compromised telo- 
merase activity and shorter telomere length were 
detected in MSC-WRN™^⁄ (fig. S7, A and B). In ad- 
dition, chromatin immunoprecipitation-quantitative 
polymerase chain reaction (ChIP-QPCR) showed 
binding of WRN to the H3K9me3-enriched cen- 
tromeric loci a-Satellite (a-Sat) and Satellite 2 
(Sat2) (Fig. 3A) (16). Depletion of WRN resulted 
in an increase in centromeric y-H2AX signal and 
a loss of H3K9me3 from o-Sat and Sat2 loci ac- 
companied by up-regulation of transcripts from 
these sequences (Fig. 3, A and B, and fig. S7C). 
Coimmunoprecipitation (Co-IP) analysis revealed 
WRN as part of a complex containing the major 
histone methyltransferase for H3K9me3—SUV39HI1, 
HP1o, and LAP2£, a nuclear envelope component 
that recruits heterochromatin via anchoring to 
HPla (Fig. 3C and fig. S7D) (17). These observations 
suggest a role for WRN, together with SUV39H1 
and HPia, in the stabilization of heterochromatin. 

We next tested whether disorganization of het- 
erochromatin could contribute to accelerated 
cellular senescence. Knockdown of SUV39H1 or 
HPlo in wild-type MSCs led to a reduction of 
overall H3K9me3 and induction of cellular 
senescence, as assayed by Western blot, SA-B-gal 
staining, and p16 expression (Fig. 3D and fig. S8, 
Ato D). On the contrary, overexpression of HPla 
up-regulated H3K9me3 levels and repressed cel- 
lular senescence in WRN-deficient MSCs (fig. 
S8, E to H). To confirm these observations, we 


Fig. 1. WRN-deficient MSCs exhibit phenotypes 
associated with premature cellular senescence. 
(A) Morphology and immunofluorescence (IF) an- 
alyses of pluripotency markers in ESCs. Scale bar, 
100 um and 10 um, respectively. (B) Western blot 
analysis of WRN expression in ESCs using anti-WRN 
N-terminal (ab200) and C-terminal (SC-5629) anti- 
bodies. (C) Growth curve analyzing the cumulative 
population doubling of MSCs. (D) Senescence- 
associated (SA)—B-gal staining in passage 1 (P1) 
and P5 MSCs. Scale bar, 50 um. (E) Quantitative 
RT-PCR analysis of the indicated genes in P1 and 
P5 MSCs. Transcript levels were normalized to MSCs- 
WRN*’* group. Genes with greater mean value are 
color coded toward red. (F) Photon flux from mus- 
cle of aNOD-SCID mouse transplanted with MSCs- 
WRN** (left) and MSCs-WRN~~ (right) expressing 
luciferase. All data are represented as mean + SEM. 
*P < 0.05, **P < 0.01, ***P <0.001 by Student's 
t test; n = 3 independent experiments. 
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Fig. 2. Epigenomic analyses of WRN-deficient MSCs. (A) Left: Chromatin 
structure of MSCs shown by Hoechst 33342 staining of the nucleus. Scale 
bar, 5 um. Right: CV value of nuclear Hoechst staining intensity used to 
evaluate the heterogeneity (pixel-to-pixel variation) of Hoechst intensity. (B) 
Immunofluorescence analyses of LAP2B expression in MSC-WRN*/* and 
MSC-WRN~“~ at P5. Arrowheads denote abnormal nuclei with decreased 
LAP28 expression (percentage of LAP2B-positive nuclei in corner). Scale bar, 
10 um. n = 3 independent experiments. (C) Enrichment of H3K9me3, 
H3K27me3, H3K4me3, and 5mC on the gene bodies and 21 kb upstream of 
TSS (transcription start site) and 21 kb downstream of TTS (transcription 
terminal site) regions in the human genome. (D) Sketch map of “H3K9me3 


mountain” distribution over 23 chromosomes. The blue lines indicate 73 
“H3K9me3 mountains” present in MSCs-WRN*’*, whereas 48 (65.8%) of 
them are localized within 5 Mb regions around the telomeres or centromeres. 
The red arrowheads indicate 28 “H3K9me3 mountains,” which are lost in 
MSCs-WRN~’~. The circles indicate the centromeres of chromosomes. (E) 
Representative images showing two “H3K9me3 mountains” on chromosome 
2 in the subtelomere or subcentromere regions in P5 MSCs-WRN~’~ and 
MSCs-WRN*7*. Two biological replicates of each sample are presented. Black 
square denotes the centromere; red rectangles denote the position of the 
presented subtelomere and subcentromere regions, respectively. All data are 
represented as mean + SEM. ***P < 0.001 by Student's t test. 
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Fig. 4. SUV39H1'"524« mutant MSCs exhibit defective nuclear envelope and heterochromatin, as 
well as phenotypes of premature cellular senescence. (A) As described in Fig. 2A (left), Hoechst 
staining images of the nucleus; right, CV value of nuclear Hoechst staining intensity used to evaluate the 
heterogeneity (pixel-to-pixel variation) of Hoechst intensity. (B) Immunofluorescence analyses of LAP2B 
expression in MSCs. Arrowheads denote the abnormal nuclei with decreased LAP28 (percentages of 
normal nuclei presented at corner). Scale bar, 20 um. n = 3 independent experiments. (C) Western blot 
analysis of the indicated proteins in MSCs. (D) SA-B-gal staining in MSCs at P5. Scale bar, 50 um. n = 3 
independent experiments. (E) Western blot analysis of the indicated proteins in human primary MSCs 
derived from old and young healthy individuals at P4 (see table S4). All data are represented as mean + 


SEM. ***P < 0.001 by Student's t test. 


generated pluripotent ESCs-SUV39HI'"""* lines 
harboring catalytically inactivated endogenous 
SUV39H1 (fig. S9, A to D). Upon differentiation, 
MSCs-SUV39H1'"*** displayed drastic nuclear 
structural and chromosomal changes, loss of 
INM proteins LAP28 and LBR, decreased levels 
of H3K9me3 and HPla, up-regulation of centro- 
meric repetitive sequence transcription, and co- 
ordinated transcriptional down-regulation of 
centromere-packaging components (Fig. 4, A to 
C, and fig. S10, A and B). MSCs-SUV39HT?™"* re- 
capitulated premature aging phenotypes observed 
in WRN-deficient MSCs, including retarded cell 
growth and accelerated cellular senescence deter- 
mined by SA-f-gal staining (Fig. 4D and fig. S10, 
C to E). High expression of SUV39H2, a germline- 
specific histone methyltransferase, and/or other 
factors may functionally compensate for SUV39H1 
deficiency in ESCs (fig. SIOF) (18, 19), where upon 
inactivation of the WRN-SUV39H1 axis, no dis- 
cernible heterochromatin change was observed 
(figs. SIOG and S8A). It should be noted that 
MSCs-SUV39H1'"** exhibit neither increased 
y-H2AX (P = 0.773) and phosphorylated ATM/ATR 
substrates (P = 0.279), nor telomere attrition (figs. 
S10H and S7, A and B). These results indicate that 
heterochromatin destabilization promotes pre- 
mature aging in MSCs. 

Finally, we asked whether heterochromatin 
disorganization could be a common hallmark for 
physiological human stem cell aging. For this pur- 
pose, we compared the levels of heterochromatin 
marks in primary dental pulp MSCs derived from 
six young (7- to 26-year-old) and six old (58- to 72- 
year-old) individuals (fig. S10I and table S4) (20). 
Amarked down-regulation of WRN protein asso- 
ciated with a decrease in H3K9me3, HPlo,, SUV39HI, 
and LAP2ß levels in MSCs derived from old 
individuals (Fig. 4E). Therefore, specific hetero- 
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chromatin changes may underlie both patholog- 
ical and physiological MSC aging. 

In summary, we have found that WRN pro- 
tein, besides its role in DNA repair, functions 
to safeguard heterochromatin stability (fig. S11). 
Our results reveal that the progressive hetero- 
chromatin disorganization observed in WRN- 
deficient MSCs underlies cellular aging, but 
more extensive studies are needed to examine 
its role during physiological aging. The meth- 
odologies and observations introduced here 
may be used and extended toward the system- 
atic study of other age-associated molecular events 
with relevance to human aging and age-related 
disorders. 
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